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ABSTRACT 

The galaxy cluster SPT-CL J0205-5829 currently has the highest spectroscopically-confirmed red- 
shift, z = 1.322, in the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. XMM-Newton 
observations measure a core-excluded temperature of Tx — 8.7l;J gkeV producing a mass estimate 
that is consistent with the Sunyaev-Zel'dovich derived mass. The combined SZ and X-ray mass esti- 
mate of A/5oo=(4.8 ± 0.8) X Mq makes it the most massive known SZ-selected galaxy cluster 
at z > 1.2 and the second most massive at z > 1. Using optical and infrared observations, we find that 
the brightest galaxies in SPT-CL J0205-5829 are already well evolved by the time the universe was 
<5 Gyr old, with stellar population ages >3 Gyr, and low rates of star formation (<0.5MQ/yr). We 
find that, despite the high redshift and mass, the existence of SPT-CL J0205-5829 is not surprising 
given a flat ACDM cosmology with Gaussian initial perturbations. The a priori chance of finding a 
cluster of similar rarity (or rarer) in a survey the size of the 2500 deg^ SPT-SZ survey is 69%. 

Subject headings: galaxies: clusters: individual (SPT-CL J0205-5829) — galaxies: formation — galax- 
ies: evolution — early universe — large-scale structure of universe 
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1. INTRODUCTION 



The South Pole Telescope (SPT;|Carlstroin et al.|20TT| 



has recently completed a survey designed to discover all 
massive galaxy clusters within a 2500 deg^ region of the 
southern sky. High redshift galaxy clusters are valuable 
as probes of the initial conditions of the universe, partic- 
ularly the distribution of matter at early epochs. Since 
galaxy clusters are the most massive collapsed systems, 
their abundance is sensitive to the properties of the early 



matter density field (e.g 


, Lucchin & Matarrese|1988 


Co- 


lafrancesco et al. 


1989 


Mortonson & Hu 2010|) and the 



nature ot inhationary models, in addition to cosmology, 
the constituent galaxies of these clusters, which have es- 
sentially co-evolving star formation histories, are useful 
for studying galaxy formation and evolution. 

The SPT-SZ survey finds clusters via t he Sunyaev- 
Zel'dovich (SZ; 
vast majority o 



Sunyaev fc Zel'dovich|[l972| ) effect. The 
: baryonic mass ol a galaxy cluster is in 
the form of diffuse, ionized gas, known as the intracluster 
medium (ICM) , unassociated with any particular galaxy. 
Photons from the cosmic microwave background (CMB) 
are Compton scattered by the free electrons in this ion- 
ized gas. The scattered photons gain energy on aver- 
age leading to a spectral distortion of the observed CMB 
known as the thermal SZ effect. The surface brightness of 
the SZ effect is independent of the distance to the cluster 
because the SZ effect depends solely on the line-of-sight 
integral of thermal pressure of the ionized gas. Therefore 
the total SZ flux is a measure of the total thermal energy 
in the gas, which is tightly correlated to the cluster mass. 
This makes SZ surveys an efficient means for finding high 
mass clusters at all redshifts (e.g., Carlstrom et al.|2002 ). 

Observations of fine scale CMB anisotropy with th e 
SPT, Planck Satellite ( [Planck Collaboration et al.|2011|), 



and Atacam a Cosmology i'elescope (AC'i'; Marriage 



et al. 20111 have recently been used to detect massive 
clusters m large surveys of the sky. The prog ress of 
the cluster s urvey by SPT is reporte d by Staniszewski 



et all ( 2009 ); Vanderlinde et al. ( 2010 ); 
mi[ )Y |Reichardt et al.| ( |2012[ ) ^FEere 



Williamson et a 
he details ot the 

survey strategy, data reduction, and cosmological anal- 
ysis are also presented. The SPT has now completed a 
survey of 2500 deg^ in the southern hemisphere in three 
millimeter-wavelength bands. The SPT-SZ survey is es- 
sentially complete for clusters with a mass of Afsoo ^ 



5 X 10" /i^ Me 



at z > 0.3 and M^qo > 3 x lO" K^M, 



70 



at z > 1 

The SPT-SZ survey has discovered several galaxy clus- 
ters that have been spectroscopically con firmed at z > 1 



( [Brodwin et al.||2010[ [Foley et al. 
are a lso discove ring high redshift c 
'e^ 



2011). Other groups 



usters thr ough X-ray 
^^^^^^ 2004; MuUis et al. 2005; St anford et al. 
Kosati et al. 2009; Henry et al. 20 10|, or mtrared " 
Stanford et al. 2005; Brodwin et al.l|2006| 
008; Muzzin et al. 2009; Wilson ct al. 



Rosati et al 

|2006 

imaging (e 
Eisenhardt etTaT 



2009,,Papovich et al..2010, .Brodwin et al.,2011,,Stanford| 
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et al.|2012 ). However, these techniques do not have the 
benefit ol simultaneously providing the nearly redshift- 
independent mass selection and wide area coverage of the 
SZ technique. 

The rich and dense environments of galaxy clusters 
can also be exploited to study galaxy evolution. They 
provide a simple way of identifying large populations of 
galaxies that have similar formation histories. As the 
redshift of the clusters in the sample increases, earlier 
phases in the process of galaxy evolution process are ob- 
served. With these observations, scenarios of galaxy for- 
mation can be tested with less temporal extrapolation. 
This also perhaps has ramifications for the hierarchical 
formation scenario, as the ages of the bright elliptical 
galaxies may be related to the epoch of the final assem- 
bly of the clu ster as suggested by s emi- analytical simula- 
tions such as 'Dubinski (1998) and jBoylan-Kolchin et al" 
( [2009| ). Recent observations of high redshitt dense envi- 
ronments point to a period of prodigious st ar formation 
( [Papovich et al.[[2012[ [Snyder et a"L[[2012[ ) at z > 1.3, 
after which the bulk of the galaxy stellar mass build up 
likely arises from dry mergers (without significant star 
formation). However, studies of lower redshift ellipti- 
cal galaxies suggest that stellar age is minimal l y affe cted 
by environment densities, e.g. Thomas et al. (2010). A 
larger data set at high redshift is probably required to 
reconcile these. 

SPT-CL J0205-5829 was first identified as a cluster in 



Reichardt et al. (2012 R12), which describes a catalog 
of 224 cluster candidates discovered in the first 720 deg^ 
of the 2500 deg^ SPT-SZ survey. SPT-CL J0205-5829 
was detected with a signal-to-noise ratio (S/N) of 10.5 
in the SPT data. Initial deep optical follow-up observa- 
tions showed no obvious overdensity of galaxies in griz 
images, but additional infrared and Spitzer photometry 
confirmed the presence of extremely red clustered galax- 
ies consistent with a redshift z > 1.3. Optical spec- 
troscopy of member galaxies confirmed that the cluster is 
at z = 1.322. X-ray observations with XMM-Aeuiton re- 
vealed a luminous and extended X-ray source. Although 
SPT-CL J0205-5829 is not the most massive SPT cluster, 
it is the highest-redshift SPT cluster that has been con- 
firmed by spectroscopy to date, and potentially the most 
massive galaxy cluster known at redshif t z > 1. 2 (pre- 
vious to this was XMM2235 from .Rosati et al.[ (|2009 1 
at z=1.39), and second most massive at z > 1 (the 
most massive being SPT-CL J2106-5844 from [Foley etal] 
( [2011 1 at z=1.13). 

We present our initial detection and follow-up observa- 
tions of SPT-CL J0205-5829 in Section [2| In Section [3| 
we show that SPT-CL J0205-5829 is a massive high re¥ 
shift galaxy cluster with a population of normal passively 
evolving galaxies. We then briefiy discuss the implica- 
tions of the existence of such a massive, evolved cluster at 
z > 1.3 in Section [1 We summarize and conclude in Sec- 
tion [5] Except where otherwise stated, we assume a flat 
ACDM cosmology with ili\i—0.3 and ho— 0.7 throughout 
this paper. M500 masses are deflned as the mass enclosed 
in a spherical region which has a density 500 times the 
critical density of the universe. At z=1.322, 1 Mpc sub- 
tends 2.0 arcminutes and the age of the universe is 4.66 
Gyr. 



SPT-CL J0205-5829 at z = 1.32 
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OBSERVATIONS, DATA REDUCTION, 
FINDINGS 



INITIAL 



2.1. Millimeter Observations by The South Pole 
Telescope 

SPT-CL J0205-5829 was initially discovered in the 
SPT-SZ survey and reported in R12, as part of the clus- 
ter catalog identified from the 720 deg^ surveyed during 
the 2008-2009 SPT observing seasons. The survey strat- 
egy and data an alysis are detailed in the previous SPT- 
SZ survey papers Staniszewski et al. j|2009 1 , iVanderlinde 
■ ■ \5wll), and R12. 



et al (20101, Wi 



Staniszewski et al. 
iiamson et al. ( |2011 

i'he SP'i'-yZ survey was completed m November 2011 
and covers an area of 2500 deg^ in three frequency bands 
at 95, 150, and 220 GHz. 

As described in R12, cluster candidates were identified 
using a multi-band m atched-filter approac h, similar to 



that first described by Melin et al. (20061. The signif- 



icance of a cluster detection (maximized across spatial 
filter scales and position in map), ^, is used to identify 
cluster candidates. For the survey field containing SPT- 
CL J0205-5829, only the 95 and 150 GHz data was used, 
the SPT maps have noise levels of 45 and 16 /iK-arcmin 
in CMB temperature units at 95 and 150 GHz, respec- 
tively. In this data, SPT-CL J0205-5829 was detected 
with ^ = 10.5 and is among the 5% most significant de- 
tections in the R12 catalog. An image of the filtered SPT 
map is shown in Figure [l] 

2.2. Optical and Infrared Imaging 

We obtained griz imaging using the M0SAIC2 imager 
on the CTIO 4 m Blanco telescope on UT 18, 25 July 
2010 and UT 4 July 2011 with mediocre to bad seeing 
(1.1 to 2.2") and occasional light clouds. Total integra- 
tion times were 300, 300, 2350, and 1050 seconds to 10-cr 
point source depths of 23.8, 23.2, 22.2, and 21.1 AB mag- 
nitudes in g, r, i, and z, respectively. We also acquired 
1800 seconds of deep i-band imaging of SPT-CL J0205- 
5829 on UT 31 January 2011 with the Ina mori Magellan 
Area! Camera and Spectrograph (IMACS; Dressier et al. 
20061 on the Baade Magellan 6.5 m telescope to 24.0 
AB magnitude depth in mediocre seeing (1.2 to 1.5"). 
The obser vation strategy an d reduction proc edure is de- 
scribe d in [High et al.| ( |2010[ HIO), Wll, and |Song et"aL 
(20121 S12) using the PHDTPIPE pipeline (Rest et al. 

'T-CL J0205-5829 was also observed with the NEW- 



FIRM imager ( Autry et aLl2003[) at the CTIO 4 m Blanco 
telescope on Ui' 6 INovember 2010. Data were obtained 
in the Kg filter under photometric conditions with a 10- 
cr point source depth of 19.1 Vega magnitudes. At each 
dither position, six frames with 10 s exposure times were 
coadded at 18 random positions providing a total ex- 
posure time of 1080 s. NEWFIRM data were reduced 
using the FATBOY pipeline, originally developed for the 
FLAMINGOS-2 instrument, and modified to work with 
NEWFIRM data in support of the Infrared Bootes Imag- 
ing Survey (A. Gonzalez, private communication). In- 
dividual processed frames are combined using SCAMP 
and SWARP (Bertin et al. 2002 ) , and p hotometry is cali- 
brated to 2MASS ( Skrutskie et al.|2006[ ). The final image 
has a FWHM of OM". 

Infrared Spitzer/IRAC imaging was obtained in 2011 
during Cycle 7 as part of a larger program to follow up 



clusters identified in the SPT survey. IRAC imaging is 
particularly important for the confirmation and study of 
high-redshift SPT clusters such as SPT-CL J0205-5829 
where the optically faint members are strongly detected 
in the infrared. The on-target observations consisted of 
8 X 100 s and 6 x 30 s dithered exposures in bands [3.6] and 
[4.5] to 10-CT depths of 20.3 and 18.8 Vega magnitudes, 
respectively. The deep [3.6] observations are sensitive 
to passively evolving cluster galaxies down to 0.1 L* at 
1 .5. T he data were reduced ex actly as in B rod win] 



et al. ( |2010[ ), following the method of |Ashby et al. (2009). 
Brieny, we correct for column pulldown and residual im- 
age effects, mosaic the individual exposures, resample to 
0'.'86 pixels (half the solid angle of the native IRAC pix- 
els), and reject cosmic rays. 

2.3. Optical Spectroscopy 

Multislit spectroscopic observations were acquired for 
SPT-CL J0205-5829 on the 6.5-meter Baade Magellan 
telescope on UT 25-26 September 2011 using the f/2 cam- 
era on the IMACS spectrograph for a total integration 
time of 11 h ours. The strategy an d procedure were as 



described in Brodwin et al. 



(2010), with the same 300 



1/mm "red" grism and WB63D(FU500 filter, but without 
the GISMO module (in order to increase throughput). 
The galaxy target selection was based on the optical 
and infrared photometry, see SectionjSj Twenty-two 30- 
minute exposures were made in excellent to moderately- 
good seeing (0.4-0.7") using one slit mask. The resolu- 
tion of the observations, as measured from the sky lines , 



was 5.2A. In a procedure identical to Ruel et ah] (2012) 
the COSMOS reduction package was used tor standard 
CCD processing, resulting in wavelength-calibrated 2D 
spectra. The ID spectra were then extracted from the 
sum of the reduced data. Spectral features were identi- 
fied by eye from inspection of the 2D and ID spectra, and 
redshifts were then obtained by using RVSAO routines. 

2.4. X-ray Observations 

A deep X-ray observation of SPT-CL J0205-5829 
was obtained by the XMM-Newton observatory (OB- 
SID:0675010101) on UT 19-20 June 2011 using the Eu- 
ropean Photon Imaging Camera (EPIC), which consists 
of two Metal-Oxide-Silicon (MOS) arrays plus one fully 
depleted p-n (PN) junction CCD array. The total inte- 
gration times were 69 ks for the MOS arrays and 65 ks 
for the PN array. The data reduction and analysis were 
performed with SAS vll.O utilizing the XMM-Newton 
Extended Source Analysis Sof tware packagep^ (XMM- 
ESAS, e.g. |Snowden et al.|2008[) . The net clean exposure 
time is 57 an d 39 ks in the MOS/PN arra ys, respectively. 
Based on the De Luca & Molendi (2004") diagnostics we 
find a ~ 30 — 40% background enhancement in the ob- 
servation due to residual quiescent soft proton contami- 
nation. The M0S2 CCD#5 was in an anomalously high 
state and we have removed it from further analysis. 

We have also excised all point sources identified in 
the source detection step. We have visually inspected 
the excision regions and made conservative adjustments 
to their size. In particular, a point source associated 
with a bright galaxy (bluer than the passively-evolving 

http : / /heasarc . gsf c .nasa. gov/docs/xnim/xinmhp_ximnesas .html 
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31.6 31.5 31.4 31.3 31.50 31.45 31.40 31.35 

R.A. <J2000) R.A. (J2000) 



Fig. 1. — (Left) The filtered SPT-SZ significance map of SPT-CL J0205-5829. The negative trough surrounding the chister is an artifact 
of the filtering of the time ordered data and maps. (Right) Color image from IMACS i, NEWFIRM Kg, Spitzer/IRAC [3.6], with SPT-SZ 
contours overlayed in white and [3.6]-[4.5] color-selected galaxies indicated in cyan. 



model) was identified in the core region of the cluster 
(a=02:05:45.4, (5=-58:28:58.3, ~ 12" west of the X-ray 
centroid) and was removed with an excision radius of 
~ 11" (see Figure [3]). 

3. RESULTS 

3.1. Cluster Member Galaxies 

From the procedure described in S12, we measure 
a redshift based on the &pitzer IRAC photometry of 
z = 1.30±0.12 (see Fi pure[2| . We fit a model of passively- 



evolved galaxies from Bruzual & Chariot ( 2003[ BC03) to 
the data to determine the redshitt. I'he optical data were 
not deep enough to offer any additional constraint to 
the redshift exce pt fo r the brightest of the cluster mem- 
bers (see section 3.3). The redshift estimator identified 
32 galaxies with iRAC [3.6]-[4.5] colors consistent with 
this redshift (a 3.5cr overdensity compared to the back- 
ground), shown in Figures [1] and [2] From this list, we 
designed a multi-slit mask for the TMA CS spectroscopic 
observations described in Section |2.3[ filling the mask 
with other targets, identified as galaxies with bluer col- 
ors relative to the model in the i-band and Spitzer data. 

3.2. Spectroscopy 

Redshifts and other spectroscopic properties of mem- 
ber galaxies are listed in Table [l] Of the 47 slits designed 
into the mask, one spectroscopic member was identified 
from an [O II] emission line, and 8 from Ca H&K. Figure 
[3]shows the spectra of the 9 cluster members. The bright- 
est cluster galaxy (BCG), defined as the brightest cluster 
member in [3.6], is at a redshift oi z — 1.3218 ± 0.0005, 
and the combined robust (biweight) redshift of 9 clus- 
ter members is 2; = 1.322^q qq2. We do not calculate 
the velocity dispersion due to the overwhelming intrin- 
sic uncertainty in the derived mass estimates with <15 



members ( Saro et ar]|2012 ). 

We also estimate the star formation rate (SFR) for 
each cluster member from the integrated [O II] flux which 
was corrected for gala ctic extincti on (reddening) using 
the dust map from Schlegel et al. (01398) and scaled to 
match the i-band magnitude from IMACS imaging. We 
do not correct for source dust extinction as we lack a 
well constrained NUV-Blue continuum measurement for 
most of these galaxies. This is also consistent with our 
derived extinction fr om S EP fits of the 4 brightest central 
galaxies (see Section 3.3). We measured the continuum- 
subtracted flux centered on the [O II] wavelength with a 
bin width of 8 A (320 km/s) and converted to luminosity 
using the cluster redshift. The SFR was estimat ed from 
the JO II] l uminosity using the scaling law from Kenni- 
The measured [O II] flux and SFR (or 3cr 
s) are given in Table fl] 

3.3. Brightest Central Calaxies 

We selected the brightest central galaxies to be the four 
brightest galaxies consistent with the [3.6]- [4.5] model, 
within 1 arcminute of the SZ center. We then use an 



cutt| ( |1998 [) 
upper limi 



analysis similar to the Rosati et al. (2009) spectral en- 
ergy distribution (SED) fitting procedure. To constrain 
the star formation history (SFH) of each of these galax- 
ies, we fit an exponential-burst stellar population SED 
model at solar metallicity and Chabrier IMF from BC03 
to the available photometry (see Figure [4]), including 
magnitude lower limits, fix the redshift at z = 1.322, 
and_add a source reddening model from |Calzetti et al.| 
(2000). From the fit parameters, we calculate the rest- 
frame K-band luminosity, stellar mass (and correspond- 
ing stellar mass-to- light ratios), and age. The uncer- 
tainties for these parameters are from the 68% con- 
fidence intervals in the multi-dimensional sampled grid 
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TABLE 1 

Spectroscopic members of SPT-CL J0205-5829 





R.A. 


Doc. 






Principal 


[Oil] flux'' 




ID 


{J2000) 


{J2000) 


z 


&z^ 


Spectral Feature 


(xlO^^^erg/cm^/s/A) 


(A^o/yr) 


J020556.50-582730.7 


02:05:56.50 


-58:27:30.7 


1.3219 


0.0007 


Ca H&K 


2.5±1.5 


0.42±0.28 


J020548.26-582848.4'= 


02:05:48.27 


-58:28:48.4 


1.3218 


0.0005 


Ca H&K 


1.4±0.9 


0.24±0.16 


J020547.70-582855.5 


02:05:47.70 


-58:28:55.6 


1.3239 


0.0003 


Ca H&K 


1.2±0.5 


0.19±0.11 


J020547.22-582901.5 


02:05:47.23 


-58:29:01.6 


1.3223 


0.0005 


Ca H&K 


<6.7 


<1.40 


J020546.57-582907.1 


02:05:46.57 


-58:29:07.1 


1.3230 


0.0005^^ 


Ca H&K 


<4.0 


<0.80 


J020543.00-582936.4 


02:05:43.00 


-58:29:36.4 


1.3119 


0.0005 


Ca H&K 


<3.4 


<0.63 


J020543.07-582956.9 


02:05:43.08 


-58:29:56.9 


1.3186 


0.0002 


Ca H&K 


<1.6 


<0.28 


J02054L06-583006.9 


02:05:41.07 


-58:30:07.0 


1.3210 


0.0007 


Ca H&K 


<5.3 


<0.90 


J020557.04-582713.9 


02:05:57.04 


-58:27:13.9 


1.3106 


0.0007 


[O II] 


14.8±4.6 


2.5±1.1 



^ Rcdshift errors arc twice those given by RVSAO. 

^ Based on integrated [Oil] flux within 4A of the line peak with no source extinction correction. Upper limits are 3cr. 
" BOG. 

^ The RVSAO error was unphysically small and was adjusted up to a typical value. 



and checked by bootstrapping this procedure hundreds of 
times and were in good agreement. These parameters are 
presented in Table [2] We find that aU models give well- 
constrained K-bandTuminosities, mainly because the ob- 
served Spitzer [4.5] filter corresponds to 2 /im in the rest 
frame. The 4000A breaks are mainly constrained by the 
deep IMACS i-band measurement. 

The rest-frame K-band luminosity of the brightest 
galaxy, at L ~ 4 x lO^^L©, is typical for BCGs in 
similar-sized clusters at z <0.25 based on previou s X- 
ray pl aarsma ct al. 20 Toj) or optica l cluster studies (Lin 
fc Mohr 2004, ,Popesso^ al.||2007[ [Brough et aT]|200gr 
an d sma!ller stuH i es ext ending to higher redshitts (z <1 ) 
by Whiley et al.| ( 2008 ) . The derived stellar mass is also 
consistent with other studies of BCGs from X -ray sam- 
ples a t similar cluster masses and redshifts (Stott et al. 
20101. The derived ages from the BC03 model hts listed 
in I'able [2] suggest that the stellar mass of these brightest 
galaxies had formed by the time of the observed epoch 
was probably complete by redshift 2 or 3, or perhaps ear- 
lier, also consistent with previous studie s of stellar ages 



of cluster galaxies at high redshift (Collins et al. 2009 



Henry et al. 2010). However this does not rule out the 
scenario of dry merging" hierarchical build up of these 
galaxies between redshift 1.3 and 3. 

3.4. iV7i? Luminosity Function 

As a further check on the cluster galaxy properties, we 
measure the observed [3.6] (roughly rest H-band) lumi- 
nosity function of galaxies with [3.6]- [4.5] colors consis- 
tent with the BC03 model from our initial redshift esti- 
mate. Galaxies selected as cluster members are within 
1 Mpc (physical distance) of the SZ-derived center and 
have [3.6]-[4.5] colors within 2a (based on each galaxy's 
photometric uncertainty) of the BC03 model. We then 
measure the number density in 0.4 magnitude bins from 
the brightest cluster galaxy to 1 magnitude brighter than 
the measured lOcr magnitude limit (to reduce any sys- 
tematic errors due to incompleteness in the catalog). 
Field galaxy contamination was corrected by measuring 
the same quantity outside of the 1 Mpc aperture and 
subtracting. We used the Schechter luminosity function. 



$(m) = 0.4;n(10)$*10"°*^("+i)ea;p(-10" 



where fi = m — m* and allowed a, and m* to vary. 
The final derived parameters and uncertainties are from 
the least squares fit to the data and bootstrapping the 
whole procedure thousands of times from the catalog se- 
lection stage. We found the [3.6] best fit parameters are 
$*=2.73 ± 0.31 arcsec-2, a=-1.02 ± 0.11 and m*=16.58 
± 0.29 (Vega), which are roughly consistent with our 
previous model assumptions of a^-1.0 and m*=17.09 at 
this redshift, calculated from the evolving stellar pop- 
ulation BC03 models, normalized to the Coma cluster 
luminosity function (sec HIO for a discussion). Recent 
measurements of the luminosity function in evol ved z > \ 
clusters fin d a similarly flat faint-end slope, a { IMancone 
et al.|[20T2 ). In contrast, less evolved high-redshitt clus- 
ters have a paucity of faint galaxies, indicated by a shal- 
lower faint-end slope (e.g. Rudnick et al. 2012 Lemaux 



et al. |2012| ) 



I'ms best fit luminosity function also corresponds to 
a richness measurement of NgQi=47±4 using the HIO 
procedure (integrating the luminosity function down to 
m*-f 1 within a 1 JVIpc physical radius of the BCG) and is 
consistent with the HIO sample of SPT-SZ clusters which 
are drawn from the same SPT-SZ significance although 
sampled at a different wavelength (observed i-band). 

3.5. SZ Mass Estimate 

We use an SZ mass estimate as described in R12 
(20111), which is calculated from the 



and [Benson et al 
rk( 



-0.4,)^ 



(1) 



Markov chain Monte Carlo (MCMC) method using avail- 
able CMB, BAG, SNe, and SPT^l (from the R12 cluster 
sample) data. The masses reported are posterior esti- 
mates based on the probability density function using the 
C and redshift for SPT-CL J0205-5829, marginalized over 
uncertainties in the SZ and X-ray (Yx) observable-mass 
scaling relations and cosmology. In Table |3j we quote 
mass estimates with and without a Bayesian prior as- 
sumption on the underlying population of clusters. The 
expected bias on the fiat-prior mass estimate is related 
to Eddington bias and affects the SPT-CL J0205-5829 
mass estimate at the ^ 10% level. This bias is due to the 
steeply falling mass function which makes it more likely 
for SPT-CL J0205-5829 to be a lower mass cluster that 
scattered up, than a higher mass cluster scattering down. 
The total uncertainty in mass (~ 20%) is dominated a 
combination of the intrinsic scatter and the uncertainty 
in the normalization of the SZ-mass scaling relation. In 
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TABLE 2 

Brightest Central Galaxy Parameters 



Galaxy 




Mk^ 






age*^ r' A^^ 




(10"Mq) 


(Vega) 


(IQiiLo) 


(Mq/Lq) 


(Gyr) (Gyr) 


A 


3.5±0.5 


-25.51±0.05 


3.3±0.2 


1.1±0.2 


4.5±0.5 O.litO.l O.OitO.l 


B 


2.9±0.5 


-25.31±0.05 


2.7±0.2 


1.1±0.2 


4.5±0.8 O.litO.l O.OitO.l 


C 


2.3±0.4 


-25.07±0.12 


2.2±0.3 


1.1±0.4 


4.5±0.6 0.1±0.2 0.1±0.2 


D 


1.9±0.3 


-25.19±0.08 


2.4±0.2 


0.8±0.2 


2.8±0.9 0.1±0.2 O.OitO.l 



^ Model Initial Stellar Mass at t— tform 
^ Absolute rest-frame K Magnitude at t— tobs 
Rest-frame K Luminosity at t— tobs 
Stellar Mass to K Luminosity ratio 

^ tobs-tform 

^ e- folding timcscalc for SFRoc exp(t/r) 

S Rcst-framc V-band extinction in magnitudes 




0.5 1.0 L5 

Red Sequence Redshifl 



"TP 



i 1 





17 

[3.6] mag 



Fig. 2.— The cluster member finding for SPT-CL J0205-5829. 
The top panel shows the excess of surface density of galaxies (above 
background) consistent within 2a of the SED model as a function of 
redshift. The RMS level of the overdensity is shown in dotted red. 
The peak overdensity is at 2: = 1.30 at 3. Scr above the background 
RMS. The middle panel shows the color-magnitude diagram of all 
objects within 2 arcminutes of the SPT center coordinates with red- 
filled black circles showing the selected galaxies from the passively- 
evolving model at z = 1.30 (black solid line) and brighter than 
m*+l. The inferred model m* is shown as a dotted vertical line 
and the [3.6] magnitude limit is shown as the dashed line. The 
bottom panel shows the Spitzer [3.6] galaxy luminosity function 
for SPT-CL J0205-5829. The dotted vertical line shows the best 
fit m* in [3.6], and the dashed line shows the [3.6] magnitude limit. 



Table |3j the mass estimates are given as M500 , defined as 
the mass within a radius in which the cluster has a den- 
sity 500 times the critical density of the universe. We 
can convert between this M500 and M200 with respect 
to Pmean, defined as the mass within a radius in which 
the cluster has a density 200 times the mean density 



of the univ erse, by assuming an NFW profile ( N avarro 



et al. 



et ai. 



1997|) and the mass-concentration relation ^ Dutty 
Using this conversion, the M200 masses are 



a factor of ^1.8 times larger, such that the unbiased SZ 
mass-estimate is M2oo=(8.7± 1.8) x lO^^'/ifg^MQ 

3.6. X-ray spectroscopy with XMM-Newton 

We estimate the X-ray physical parameters of SPT- 
CL J0205-5829 using an iterative process over the cluster 
radius. We measure the core-excised X-ray temperature, 
Tx, within rsoo, defined as the radius inside which the 
mass density is higher than 500 times the critical density 
of the universe. We iterate over values of r^QQ so that 
the measured Ty maintains consiste ncy with the M — T 
relation from iVikhlinin et al.| ( [20091 V09). 

For each value of rsoo, we extract spectra and redistri- 
bution and ancillary response files. We excise all detected 
point sources from both the source and background re- 
gions as well as the central r < 0.15 rsoo cluster core re- 
gion. Give n th e significant residual quiescent contamina- 
tion (Sect. |2!4| , we opt to use a local background model 
in the fitting procedure. For each camera, we subtract 
a background spectrum extracted from an annulus cen- 
tered on the cluster between 160" and 320" in radius. 
These radii were selected based on the cumulative count 
rate profiles so that the annulus is not contaminated 
by cluster emission while still lying on the same IVIOS 
chips as the source. The total number of background- 
subtracted source counts is ~ 5500 for all three cameras. 
We use Xspec vl2.5 to fit the spectra with a MeKaL 



model (Ivlewe et al. 1985 Kaastra et al. 1992 Liedahll 



et al.|l9U^ I using C-statistics on minimally bmned spec 



tra (i e., binning only channels to obtain > 1 counts/bin). 
From this fit to the spectrum, we measure the X-ray tem- 
perature. 

We then use the measured X-ray temperature and the 
redshift from optical spectroscopy to infer a M500 mass 
from the V09 M — T relation, which we also convert 
to a corresponding rsoo value. Given this new value of 
^500 J we iterate on this process until two successive rsoo 
estimates differ by < 2.5" (equal to the bin size of our X- 



SPT-CL J0205-5829 at z = 1.32 




Rest-frame Wavelength (Angstroms) 
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Fig. 3. — (Left) Optical (r/i) and Spitzer/IRAC [3.6] image showing the galaxies confirmed by spectroscopy and overlayed by the XMM- 
Newton X-ray contours. Cyan shows the [O II] detection shown to the right (bottom panel), and red shows the members identified with Ca 
H&K features shown to the right. The white circle shows the X-ray point source position. The frame subtends 4.5x3.4 arcminutes. (Right) 
Spectra of all 9 member galaxies with the [O II] and Ca H&K features indicated in vertical dotted lines. Despite the long exposure time 
on Magellan, these features arc faint due to the red color of passively galaxies. 



ray images). This criterion was reached in four iterations 
and we have verified that the final solution is independent 
of the initial rsoo value. The rsoo radius is 710 kpc(~ 85") 
and the excised core region (i.e., O.lSrsoo) has a radius 
of ~ 13" (roughly twice the PSF FWHM). 

The final spectrum is displayed in Figure [5] The best 
fit temperature is Tx =8.7;to g keV. This corresponds to 
a mass of Af5oo=(5.2± 1.3) x IO^^HyqMq using the rela- 
tion from V09. There is no evidence that SPT-CL J0205- 
5829 is unrelaxed from the X-ray morphology or the 
galaxy distributions, either on the sky or in the veloc- 
ities, therefore we use the standard M — T relation from 
V09 and make no corrections based on the dynamical 
state of the cluster. In the Tx-based mass estimate, we 
include the statistical uncertainty in the measurement of 
Tx, uncertainties in cosmology, and assume a 20% in- 
trinsic scatter in the M — T relation, as noted by V09. 

We constrain the mean [Fe] abundance to 
Z=O.26±O.15Z0 which is consistent with a typ- 
ical mean metallicity for the ICIVI in a mas- 
sive galaxy cluster (Z^ 0.3Zf:i) at lower redshifts 



iVlaugnan et al. 



IVlatsumoto et al. 



^oom 



[Tozzi et al. 





2003 


m 


''soo 



is Lx(0.5-2.0 keVy^ (3.91 ± 0.05) x 10'*'' erg s'^ in 
the rest frame. 

We also note that we detect the Fe K line clearly in 
PN and JVI0S2 and more weakly in JVIOSl. If we allow 
the redshift to vary during fitting the best fit value from 



TABLE 3 

Mass Estimates for SPT-CL J0205-5829 



Observable 


Measurement 


M5oo(1O"'ito1^0) 


SZ 5 


10.5 


4.8 ± 1.0 


SZ S. (fiat prior) 


10.5 


5.2± l.l 


Tx 


8.7+1° keV 


5.2 ± 1.3 


Combined 




4.8 ±0.8 



Note. — Uiiboldcd masses indicate Eddington-biased 
mass estimates, calculated using fiat priors on mass. Note 
that the SZ (untargcted) flat-prior measurement suffers 
from a considerably different Eddington bias than the 
(targeted) fiat-prior estimate. Tlic Tx mass estimate and 
the unbiased SZ ^ mass estimate were used to generate the 
combined mass estim ate based on their probability distri- 
butions, (see Section |3.7[ . 

the joint fit is z=1.39 ± 0.02, which is ^ 5% larger than 
the optical spectroscopic redshift of z = 1.322^Q'gQ2- The 
redshifts derived from individual cameras are Ifi consis- 
tent with the joint PN+M0S1-I-1VI0S2 fit (except IMOSl 
which gives a slightly lower redshift). This is one of 
the highest redshifts me asured from X-ray spectra (cf., 
Lloyd-Davies et al.||2011[ ). 



3.7. Combined Mass Estimate 



We follow Foley et al. (2011 ) and calculate a joint es- 
timate using the SZ mass and X-ray mass estimates of 
SPT-CL J0205-5829. We assume the uncertainties are 
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Fig. 4. — Best BC03 models for the four brightest galaxies in the central region of SPT-CL J0205-5829. The left panel shows thumbnails 
of 6 observed filters (r, i, z, Ka, [3.6], and [4.5]) in the central ~30" from the brightest galaxy (A), with the four galaxies labeled in the Kg 
image. The right panels show the measured photometry in red with the best fit BC03 model for each overplotted. 
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Fig. 5.— The X-ray spectrum of SPT-CL J0205-5829 (black: 
PN, green: MOSl, red: MOS2). 



uncorrelated between the two masses. This allows for a 
more straightforward evaluation of the posterior proba- 
bility distribution function (PDF), 



P{M\C,Tx) (X P{M)P{C\M)P{Tx\M), 



(2) 



et al.||2008] ), P{£,\M) is the flat-prior SZ mass estimate 
FDi"', and P{Tx\M) is the flat-prior Tx mass estimate 
PDF. As calculated in Sections 3.5 and 3.6, we use the 
X-ray and SZ mass estimates derived from the observ- 
ables Tx and ^, respectively, which were marginalized 
over uncertainties in their scaling relations and cosmol- 
ogy. We find a combined, unbiased, mass estimate to be 
M5oo={4.8 ± 0.8) X lO"/ifo^M0. Converting to M200 as 
above gives M2oo=(8.8 ± 1.4) x lO"' 



4. DISCUSSION 

The galaxy members in SPT-CL J0205-5829 were iden- 
tified via [3.6]-[4.5] color, and the significant overdensity 
of these is how this cluster was initially confirmed after it 
was identified by the SZ effect. The measured overden- 
sity (richness) of these galaxies is consistent with other 
SPT-SZ clusters, a sample that has a median mass of 
M500 ~ 3.3 X lO^/ifo^ Mq. From SED fitting, the BCG 
and 3 other bright central galaxies have luminosities and 
stellar masses typical of central galaxies in clusters of 
similar mass at lower redshift and have derived stellar 
population ages greater than ~3 Gyr. This suggests that 
most of the eventual stellar mass in these galaxies are al- 
ready present at z —1.3 and the vast majority of these 
stars were formed by z ^3. The actual assembly scenario 
of these galaxies cannot be constrained yet. 

The quiescent SEDs and the amount of [O II] in the 
spectra of the central galaxies suggest there is very lit- 
tle ongoing star formation (<O.5M0/yr) in the center of 
SPT-CL J0205-5829, meaning there is no strong cooling 



Hu et al. 


19851 


19991 Hate 


h et a 



H eckm an et 
W7t 



al 



1989 



McDonald et al 



Cr awfor d et 



al. 
How- 



where P[M) is the Tinker halo mass function (Tinker ever, it should be noted that only a very strong cool- 



SPT-CL J0205-5829 at z = 1.32 
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ing core wou ld be discernible in [O II] ( McDonaId]|2011 
Santos et al. 2011) and a better indicator wou ld b e tig- 
'I'he central A-ray point source noted in Section 2.4 could 
possibly be a central AGN that is suppressing the star 
formation in the cluster. Alternatively, given the age 
of the cluster and the X-ray cooling time, there might 
not have been sufhcient time for a strong cooling flow to 
form. Either interpretation is consistent with the general 
lack of star formation and strong cooling flows at z > 0.5 
found in previous studi es (Santos ct al. 2008', Vikhlinin 
et al.||2007[ [Samuele et al. 2011; McDonald 20lT|), where 
the number density ol strong cooling hows increases dra- 
matically at z <0.5 while AGN and merger activity de- 
creases. We do not see any indication of a major merger 
either in the X-ray morphology or galaxy distribution. 

From the X-ray spectrum, we found that the cluster 
gas has a metallicity consistent, albeit with significant 
uncertainty, with massive clusters at lower r edshift. This 
is als o consistent with several studies (see Baldi et al. 
2012 for a review) of ICM metal abundances over a range 
of redshifts that found little or no evidence of evolution 
from z <1.4. Several studies have found that this enrich- 



ment can happe n over a timescale of 1 Gyr (Pipino & 
Matteucci|2004 ) and settle into the central region withm 
a cluster crossing time (1 Gyr). The best-fit metallic- 
ity of SPT-CL J0205-5829 would suggest that the bulk 
of the metal production could have been completed by 
z ~2.5, however more X-ray observations are needed to 
say this with high statistical significance. 

The optical and infrared data have shown that the 
stellar populations of the most massive central galaxies 
are already well-evolved, suggesting that the assembly of 
these galaxies happened within the preceding 2-4 Gyr. 
It has been suggested that this timeline may depend on 
the mass of the cluster, as there is some evidence that 
the BCGs at z > 1 in lower mass clusters have not fully 
assembled (Stott et al.|2010). It may become possible to 
see this change over the full SPT mass range and may 
probe different regimes where other feedback modes dom- 
inate. Such a study of a large SZ-selected sample has the 
potential to directly measure the build up of the stellar 
mass as a function of redshift and cluster mass. 



4.1. Rarity 

Although SPT-CL J0205-5829 was included in the 
sample used in R12 for cosmological analysis, we did not 
assign a goodness of fit to the model, so it is interesting 
to quantify the probability of having found this cluster 
in the full 2500 deg2 SPT-SZ survey. We use the fuU 
2500 deg^ SPT-SZ survey area in order to avoid a poste- 
riori selection of the area in which SPT-CL J0205-5829 
was found, which could artific i ally b oost the apparent 
rarity. We follow |Foley et al. (2011) and compute the 
probability of finding a cluster at higher mass and higher 
redshift than SPT-CL J0205-5829. We do so by sam- 
pling the cosmological and scaling relation constraints of 
the CMB+BBN-l-BAO-l-HST-t-SN+SPTcL chain from 
R12 and producing a posterior statistical mass estimate 
P{M\^, z) at each step in the chain. We then com- 
pute the expected number of clusters at higher mass and 
higher redshift 



->z>M 



dN 
dM'dz' 



M' 



P{M"\^,z)dM"dz'dM', 
(3) 

whe re is the mass function as calculated follow- 

ing ^' 



Tinker et al. 



mological and sea 



(20081. The median point in cos- 
ing relation parameter space predicts 
x>z>M = 0.07 clusters at higher mass and higher redshift 
than SPT-CL J0205-582 9 in 2500 de g^ 



How e ver, as noted by 



HotchkissI (|2011j); |Hoyle et al , 

(20121; Waizmann et al. ( |2012a|b[ ), this statistic has a 
small expectation value due to the fact that it requires 
a cluster of simultaneously higher mass and higher red- 
shift than a particular object. This statistic does not 
consider the fact that many similarly rare clusters could 
exist with a slightly higher mass and lower redshift or 
lower mass a nd higher r edshif t. Instead, we follow the 
treatment of Hotchkiss (20111 and compute the proba- 
bility of finding the particular value of 5;>z>m = 0.07, 
corresponding to SPT-CL J0205-5829 for an ensemble 
of simulated 2500 deg^ surveys. We then create a nor- 
malized histogram of the resulting values of Xyz>M for 
the rarest cluster in each catalog and integrate the area 
under the curve from to the value of i>z> m for the par- 
ticular cluster in question. This statistic, unlike Xyz>M 
itself, has an expectation value of 0.5. We note that 
it depends only very weakly on the details of the sim- 
ulation or the point in cosmological or scaling relation 
space at which the simulations are performed. This met- 
ric suggests that this cluster is not at all surprising with 
a probability of 0.69 of finding at least one cluster as rare 
as SPT-CL J0205-5829 in 2500 deg^ . 

As a comparison to other more rare clusters in the 
SPT-SZ survey, using the same statistic we find 0.21 for 
SPT-CL J2106-5844, which was considered in |Foley et ah] 
(201 1| ) and 0.05 for SPT-CL J0102-4915 (ACT-CL 0102- 
4915), currently the rarest cluster in 2500 deg^ SPT-SZ 
survey. 

5. SUMMARY AND CONCLUSIONS 

We report the massive galaxy cluster SPT-CL J0205- 
5829 at z = 1.322 discovered in the first 720 deg^ of the 
SPT-SZ survey and present results of follow up observa- 
tions at optical, infrared, and X-ray wavelengths. The 
galaxy population of this cluster shows a strong red se- 
quence with a luminosity function consistent with that 
of lower redshift SZ-selected clusters. Galaxy SED fits 
to an exponentially decaying SFR stellar population, the 
[Fe] abundance from the X-ray spectrum, and the lack of 
[O II] emission in most of the optical galaxy spectra sug- 
gests that the bulk of the star formation happened at an 
earlier epoch {z > 2.5). Optical spectroscopy of 9 galax- 



ies confirms the cluster redshift at z = 1.322; 



-0.001 
-0.002 



, also 



roughly consistent with X-ray spectroscopy which gives 
z=1.39 ± 0.02. This estabhshes SPT-CL J0205-5829 as 
the highest-redshift SZ-selected galaxy cluster verified by 
spectroscopy, and the second most massive SZ-selected 
cluster known at z > 1. Based on the X-ray temperature, 
SPT-CL J0205-5829 is consistent with being more mas- 
sive than XMM2235 at z=1.39 with Tx =8.6lH from 



iRosati et al. (2009) (but the uncertainties in both tem- 
peratures are much larger than the measured difference). 
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The measured mass observables (from the SZ and X- 
ray temperature) are consistent and give a combined 
mass estimate of M500 = (4.8 ± 0.8) x lO^/ifo^Mo. Al- 
though not the most massive SZ-discovered cluster, it 
demonstrates that a cluster of this mass has enough time 
to form during the first 5 Gyr of the universe, and the 
existence of this rare object appears to be fully consistent 
with general expectations for a flat ACDM cosmological 
model. 

In general, we find that SPT-CL J0205-5829 has prop- 
erties similar to clusters with the same mass at lower 
redshift. This is extremely important in the context of 
an ultimate goal of an unbiased and low scatter mass- 
calibration of clusters over a wide range of redshifts for 
cosmological studies, and provides new insight to the as- 
sembly of the rarest and most massive structures in the 
universe. 
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